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Abstract: The prospects for a determination of the strong coupling constant as 
via scaling violations of fragmentation functions in deeply inelastic scattering are 
studied. The statistical error in the case of an integrated luminosity of 250 pb~^, 
and the theoretical errors due to the various parton density parametrizations and 
to the factorization scale dependence are estimated. 



1 Introduction 



The strong coupling constant as has been measured at HERA by means of the (2+1) jet rate ||l|. 
This particular process has the advantage that as is, up to higher order corrections, directly 
proportional to the ratio of the measured cross-sections. Another route to a determination of 
as is given by scaling violations of phenomenological distribution functions. Perturbative QCD 
predicts the scale evolution of these quantities by means of renormalization group equations 
[0,^. As a consequence, again up to higher order corrections and resummation effects, the ex- 
perimentally measurable quantities (the structure functions, and the fragmentation functions, 
which depend on factorization scale and scheme) are (symbolically) of the form A + i? In /i^. 
Here /i is the factorization scale, which is to be identified, in deeply inelastic scattering, with 
a scale of the order of the photon virtuality Q, for lack of other hard scales related to the 
leading-order process. The distribution functions contain a /i-independent term A, and the 
as dependence is only logarithmic in the factorization scale. The scaling violations are there- 
fore expected to be small, and will require large luminosity to be statistically significant. An 
as determination via scaling violations has the advantage that, in principle, no explicit model 
assumptions such as specific fragmentation models go into the measurement. In the case of 
scaling violations of structure functions^ a completely inclusive quantity is measured, and the 
theoretical basis, namely the operator product expansion, is very transparent and can be de- 
rived rigorously from light cone dominance. In the case of one-particle-inclusive processes, 
where the operator product expansion is not available, the factorization theorem of perturba- 
tive QCD (see, for example, Ref. |^ and references therein) allows the separation of the hard 
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scattering process from the non-perturbative fragmentation process. The one-particle-inclusive 
cross-section is a convolution of a mass-factorized parton-level scattering cross-section, a parton 
density and a fragmentation function: o = cxhard ® / ® -D. A possible strategy for an mea- 
surement at HERA is to perform a combined multiparameter fit of fragmentation functions 
and of the strong coupling constant to the Xi^-distribution p{xf) = (da/dxp) /o"tot (or to any 
other distribution sensitive to the fragmentation functions) of charged hadrons at two different 
scales Q. Here the variable xp is defined to be 2hL/W, where hi is the longitudinal momentum 
fraction of the observed charged hadron in the direction of the virtual photon in the hadronic 
centre-of-mass frameQ and W is the total hadronic final state energy. In leading order, xp is 
the momentum fraction of the final state current quark carried by the observed hadron. The 
total cross-section is denoted by cxtot- The strong coupling constant enters the expression for 
p in three places: (a) as an expansion parameter in the next-to-leading order expression for 
(Thard; ciiid lu the reuormalizatiou group equations of (b) the fragmentation functions and (c) 
the parton densities. The parton densities are an input to the analysis. Since they are obtained 
by a global fit, where a specific value of as is used, it is necessary to include this dependence 
as well as the variation due to the different parametrizations into the systematic errorQ. The 
next-to-leading-order one-particle-inclusive cross-section has been calculated in Ref. 0. For 
our study, we use a recent recalculation and numerical implementation described in Ref. ||^. A 
comparison of the theoretical x^-distribution with experimental data from the HI and ZEUS 
Collaborations P] has been done in Ref. . It turns out that the theoretical description of the 
experimental data is quite satisfactory. The next-to-leading order result is always within one 
standard deviation of the experimental data points except for those at very large xp, where the 
currently available fragmentation function parametrizations are not well constrained by e'^e~ 
data. 

In the next section, we describe the estimate of the various errors^ of the value of as- We 
also discuss the sensitivity to the strong coupling constant of various ranges in xp- The paper 
closes with a short summary and conclusions. 

2 Error Estimates 

To get a quantitative estimate of the dependence of the scale evolution of fragmentation func- 
tions on the employed value of as, we fix the fragmentation functions at a scale of /io = 2 GeV 
as the leading-order parametrization of Ref. |jTO|. We then evolve this input with two different 
values for Aq^j^ of 0.1 GeV {a} and 0.2 GeV {b}. The corresponding ajij'-distributions p^""^ and 
p^^^ are determined for these two sets of fragmentation functions. We now assume that the 
XiT-distributions p are measured in two different bins i, j of the factorization scale Q. The 
ratios A'"^*-'^' = p^"^^ / p^^^ for an arbitrary coupling constant as (taken at the mass of the Z boson) 
are expanded in a power series in a^, where only the linear term is kept: 



^It might be possible to reduce the dependence on the parton densities and on the not yet weU understood 
physics of the forward direction by performing an analysis in the Breit frame. I thank N. Brook and T. Doyle 
for remarks concerning this issue. 

■^Alternatively, parton density parametrizations with varying values of the strong coupling constant |^ can 
be employed. 

We do not consider experimental systematic errors. 
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Based on this formula, a quantitative estimate of the statistical error of is possible. Moreover, 
by varying A^*-'^, for example by using various parton density distributions or by modifying the 
factorization scale, the impact of systematic effects can be studied. 
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Figure 1: Statistical errors of as for the individual xp-bins, for the Q-bins {!}, {2} 
{2}, {3} [---]. 



and 



To be more specific, we assume the phase space cuts and a;i7'-bins of the ZEUS analysis, 
except for the range in the photon virtuality Q, where we consider three bins: [3.16, 12.6] GeV 
{1}, [12.6, 100] GeV {2} and [100, 150] GeV {3}. To obtain explicit numerical values, we use 
the CTEQ 3L parametrization [|TT| for the parton densities (for simplicity, we work in leading 
order). The integrated luminosity is assumed to be 250 pb~^. Under the assumption of Gaussian 
statistical errors we arrive at a statistical error of as{M^) of Aaf^^ = ±0.0007 for an analysis 
based on bins {1} and {2} and of ±0.027 for an analysis based on bins {2} and {3}. The 
individual errors of for the various Xi?-bins are shown in Fig. |l|. The large error around 
xp ^ 0.08 for the large-Q bins comes from the fact that the evolution of the fragmentation 
functions around this value of the momentum fraction is quite small (for smaller values, the 
fragmentation functions increase with increasing factorization scale, and for larger values, they 
decrease). The sensitivity of the cross-section to a variation of is largest at large xp, but this 
region also suffers from small statistics of the data sample. It turns out that the full Xp-Tsnige 
is about equally important. 

As briefly mentioned already in the introduction, an input parton density has to be chosen. 
To estimate the size of this effect, we determine the spread of the results for as(M|) depending 
on the next-to-leading-order parton densities from Refs. |Tl],|12[. For the bins {1} and {2}, the 
spread is Aa^^^ = ±0.017, and for the bins {2} and {3}, it is ±0.005. Future global fits of 
parton densities including improved HERA data should reduce this systematic uncertainty. 

Perturbative QCD allows for some freedom in the choice of the factorization scale /i of the 
fragmentation functions D{z,^'^). This brings out the inherent uncertainty in the theoretical 
prediction, and can be interpreted as an effect of unknown higher order contributions. To obtain 
an estimate of this uncertainty, the ratios A are determined for the three choices Q and 
2Q of this scale. The change of cross-section has for consequence a variation in the extracted 
a,(M|) value of Aaf^^" = ±0.013 and ±0.011 for the combinations of the bins {1}, {2} and 
{2}, {3}, respectively. 
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Figure 2: The running strong coupling constant, as recently observed in high energy collider 
experiments. 

3 Summary and Conclusions 

We have studied the prospects of a measurement of the strong couphng constant in deeply 
inelastic scattering at HERA by means of scaling violations of fragmentation functions. The 
combinations of the obtained values for Aaf^*, Aa^°^ and Aa^'^^'*' are large compared with the 
present error Aa^ = 0.006 of the world average. It is therefore likely that a measurement of this 
kind will not be competitive, concerning the size of the error. Nevertheless, it is worth doing 
as an independent quantitative test of QCD and, more important, because it complements 
the other (potential) HERA measurements based on (2+1) jet rates and scaling violations of 
structure functions. The prospects for the observation of the running of as (Fig. |^) should also 
be studied in some detail. 
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